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Abstract 

It is shown that the averaged velocity profiles predicted by statistical mechanics of point vortices 
and statistical mechanics of vortex lines are practically indistinguishable for a shear flow between 
two parallel walls. 
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I. INTRODUCTION 



It is known that statistical mechanics of point vortices describes surprisingly well the 



averaged velocity profi 
chanics of vortex lines 



es of self-similar mixing layer. After development of statistical me- 



, where a vortex does not remain straight and is allowed 
to take wavy shapes in the course of motion, it appeared a concern that the above-mentioned 
feature of statistical mechanics of point vortices can be lost in three-dimensional theory. In 
this paper we show that this is not the case: the averaged velocity profiles predicted by sta- 
tistical mechanics of point vortices and statistical mechanics of vortex lines are practically 
indistinguishable for a shear flow between two parallel walls. 



II. AVERAGED EQUATIONS FOR FLOW BETWEEN TWO PLATES 

Consider a flow of ideal incompressible fluid between two parallel walls, y is the coordinate 
normal to the walls, —h < y < h, 2h - the distance between the walls. The flow is modeled 
by motion of a large number of vortices (Fig.l). Flow is periodic in x-direction and z- 
direction. In the limit of infinite period in z-direction, the averaged velocity is parallel to 
the walls, does not depend on x and has the only non-zero component, u = u(y). Assuming 
that vortices have the same intensity and the total discharge is zero, one obtains for the 
stream function of the averaged flow, ip = ip(y) (u = difj/dy), the equation 



= ±U, (2.1) 

±h 

where f(y) is the probability to find a vortex at the point y, a is the total vorticity per unit 
length in x-direction, a = —2U. In the case of point vortices (Fig. la), {3,0] 

e —P<nji(y) 

f(y) = - h :— ; (2-2) 

and equations (1), (2) form a closed system of equations. This system can be solved an- 
alytically. Parameter (3 has the meaning of inverse temperature of vortex motion. It is 
determined by the initial energy of turbulent flow. 

In the case of deformingvortex lines (Fig. lb), f(y) is expressed through the solutions of the 
eigenvalue problem, {4, 3] 
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FIG. 1: Two-dimensional flow of point vortices (a) and three-dimensional flow of vortex lines (b) 
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FIG. 2: Stream function (a) and velocity (b) profiles for fop = — 1.0 and @2D = —1-1374 



Alp — fiaipip = —\(f, 



d 
dy 



0. 



(2.3) 



±h 



A being the minimum eigenvalue. Similarly to quantum mechanics, f(y) is proportional to 
the squared solution of the eigenvalue problem: 



L h ( p 2 (y)dy' 



(2.4) 



We aim to compare solutions of (1),(2) and (1), (3), (4). 
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III. RESULTS AND DISCUSSION 



There is no reason to expect that the velocity profiles found from the two quite different 
system of equations, (l)-(2) and (1),(3),(4), coincide. Nevertheless, this turns out to be 
the case: the velocity profiles are practically indistinguishable. More precisely: for each (3 
from "3-D problem" (1),(3),(4) there is (3 from "2D problem" (l)-(2) for which the velocity 
profiles practically coincide. 

System of equations (1), (2) admits an analytical solution: 

V = -^lncosh(C 2/ ), (3.1) 

where the constant C is determined from the boundary condition, 

C 

U + — tanh(CTi) = 0. (3.2) 
For h = oo, C — —U 2 j3 and (5) becomes: 

V> = -^lncosh(f/ 2 /fy). (3.3) 

All parameters of the flow can be normalized with respect to the wall velocity: ip — > ifr/U, 
a — > a/U and (3 — > (3U 2 . The dimension of the normalized (3 is length -1 . The characteristic 
length, \j3\~ 1 , is the width of the mixing layer. 

Non-linear eigenvalue problem (1), (3), (4) was solved numerically using an iteration pro- 
cedure. At each iteration step the system is being reduced to the Sturm-Liouville problem, 

ml 

that was solved with the Priifer method (a shooting method based on oscillation [ll, l9(). For 
a given (3 of vortex line flow, (3^d, one can choose f3 of point vortex flow, (3 2 d-, in such a 
way, that the velocity profiles are practically identical. This can be seen, for example, from 
Fig. 2, where the velocity profiles and the stream function are shown for (3^b = —1.0000 and 
j3 2 D = —1.1374. The values of (3 2 d and j3 3 D for which the average velocity profiles coincide 
form a curve in the plane {^d-i^d} shown in Fig. 3. This curve was found in the following 
way. For each f3%D we seek (3 2 d by the minimization of the sum 

N 

$>h 2D -^l, (3-4) 
i=i 

where u 2 d and u^c are velocities in 2D and 3D problems, respectively, N is the number 
of the mesh points. The weight coefficients, Wi, depend on the density of the mesh and 
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FIG. 3: Two-dimensional inverse temperature @2D vs - three-dimensional inverse temperature (5^d 

the gradient of the velocity profile u. The distance between the walls, h, was chosen large 
enough for the solution (7) to be applicable. 

The parameters fiio an d (3%d have the meaning of the inverse temperature of two- 
dimensional and three-dimensional motions. The fact, that the same velocity profile 
corresponds to different /?2d and (3zd, indicates that the corresponding temperatures of 
two-dimensional and three-dimensional motions are different. The temperature of two- 
dimensional motion has a simple physical meaning: this is an average area bounded by 
the vortex trajectory. The area has orientation. Accordingly, temperature may have both 
signs. Negative temperature of the flow considered corresponds to clockwise pass of the curls 
of the vortex trajectories [2j. The physical meaning of temperature of three-dimensional vor- 
tex line motion is not known, therefore a physical interpretation of the graph /32d{@3d) is 
yet to be established. The result obtained seems an indication that for the shear flow be- 
tween parallel walls three-dimensionality does not play an important role in formation of 
the averaged velocity profiles. 
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